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Abstract. Given two measured spaces {X,fi) and {Y,w), and a third space 
Z, given two functions u{x,z) and v(x,z), we study the problem of finding 
two maps s : X — > Z and t : y — > Z such that the images s(fi) and t(i') 
coincide, and the integral J-^ u{x, s(x))dfi + Jy v{y,t{y))di' is maximal. We 
give condition on u and v for which there is a unique solution. 



Suppose we are given three goods, X, Y, and Z . They are not homogeneous, but 
come in different qualities, x G X , y ^ Y and z G Z. Goods X and Y are used for 
the sole purpose of producing good Z, which we are interested in. To obtain one 
piece of good Z, one has to assemble one piece of good X and one piece of good 
Y. More precisely, one can obtain a piece of quality z by assembling one piece of 
quality x and one piece of quality y, yielding a benefit of u (x, z) + w (y, z). Given 
the distributions and v of goods X and Y, one wishes to minimize the total 
benefit of production. 

This translates into the following optimization problem: find maps s : X Z 
and t : Y ^ Z such that s (/z) = t (v) (this is the matching condition) and the 
integral 



is maximized. 

The origins of that problem lie in the economic theory of hedonic pricing (see |S] 
for an overview). The economic aspects will be developed in another paper P]. 

Mathematically speaking, this is related to the classical optimal transportation 
problem (see the monographs 8 and 9 for accounts of the theory). Recall that 
this problem consists in minimizing the integral 



among all maps s : X Y such that s (fi) = v. Here the measured spaces {X, fi) 
and (Y, are given, as well as the function u : X x Y s- R. A seminal result by 
Brenier ^ states that, if X and Y are bounded open subsets of i?" endowed with 
the Lebesgue measure, with X connected, and u {x, y) = ||a; — y||, then there is a 
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unique solution s to the optimal transportation problem, and s is almost everywhere 
equal to the gradient of a convex function. 

Kantorovitch 7 introduced into the optimal transportation problem a duality 
method which will be crucial to our proof. Instead of proving directly existence 
and uniqueness in the optimal matching problem, we solve in section 3 another 
optimization problem, and we will show in section 4 that it yields the solution to 
the original one. This correspondence relies heavily on an extension of the classical 
duality results in convex analysis (see U). This extension has been can be found 
in 121 and [?]; for the reader's convenience, we will give the main results in section 
2. Finally, in section 5, we will give some consequences of the main result. 

From now on, X C i?"^ , Y C i?"^ , and Z C i?"^ will be compact subsets. We are 
given measures fi on X and on y, which are absolutely continuous with respect 
to the Lebesgue measure, and satisfy: 

fi {X) ^iy{Y)<oo 

We are also given functions u : fii — > i? and v : Q2 ^ R, where fli is a 
neighbourhood oi X x Z and ^1-2 is a neigbourhood oiYxZ. It is assumed that u 
and V are continous with respect to both variables, differentiable with respect to x 
and y, and that the partial derivatives D^u and DyV are continous with respect to 
both variables, and injective with respect to z : 

(1.3) Vx e X, DxU [x, Zl) = DxU [x, Z2) =^ Zl = Z2 

(1.4) \fy e Y, DyV {y, zi) = DyV (y, Z2) =^ zi = Z2 

The latter condition is a generalization of the classical Spence-Mirrlees condition 
in the economics of assymmetric information (see |2|). It is satisfied for u {x, z) = 
— provided a ^ and a ^ 1. 

Theorem 1. Under the above assumptions, there exists a pair of Borelian maps 
(s, t) with s (fj.) = t (v), such that for every (s, t) satisfying s (p.) — t (i/), we have: 

/ u {x , s {x)) dfi — / v{y,t{y))di^ < / u{x,s [x)) d^i — I v{y,t{y))dv < 00 
Jx Jy Jx Jy 

This solution is unique, up to equality almost eveywhere, and it is described as 
follows: there is some Lipschitz continuous function p : Z R and some negligible 
subsets Xq C X and Yq C Y such that, for every x ^ Xq and every y ^Yq: 

(1.5) yz^s{x), u {x, z) — p (z) < u {x,s (x)) — p {s {x)) 

(1.6) Vz^t(y), v{y,z)-p{z) > v(x,t{y)) ~p{t{y)) 

(1.7) pis differentiable at s (x) and f (y) 

If in addition u and v are differentiable with respect to z, we get, from the 
minimization (|1.5|1 and the maximization (|1.6() : 

D^u (x, s{x)) ^ D^p (s (x)) 
D,v{y,t{y))^D,p{t{y)) 

Set s(/i) —t{i') = A. It follows from the above that, for A-almost every z £ Z, 
there is some x ^ Xq and y ^ Yq such that z — s (x) = t{y), and for every such 
(x, y) we have: 

DzU (x, z) = Dzp {z) = DzV (y, z) 
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Note that there is no reason why A should be absolutely continuous with respect 
to the Lebesgue measure. 

The proof ot theorem^s deferred to section 4. Meanwhile, let us notice that we 
have slightly changed the formulation of the optimal matching problem: by setting 
w = —V we recover the original one. This change will simplify future notations. 

2. Fundamentals of u-convex analysis. 
In this section, we basically follow Carlier 

2.1. u-convex functions. We will be dealing with function taking values in RU j+cxo}. 
A function / : X — > MU {+cxd} will be called u-convex iff there exists a non-empty 

subset A C Z X M. such that: 

(2.1) \/xeX, sup {u{x,z) + a} 

(z,q)gA 

A function p : Z MU {+cxd} will be called u-convex iff there exists a non-empty 
subset B C X X R such that: 

(2.2) p{z)= sup {u{x,z) + b} 

(x.b)eB 

2.2. Subconjugates. Let f : X ^ KU{+oo}, not identically {-|-cxo}, be given. 
We define its subconjugate p : Z ^ MU {+00} by: 

(2.3) /«(z)=sup{u(a;,z)-/(a;)} 

X 

It follows from the definitions that /" is a u-convex function on Z (it might be 
identically {+cx)}). 

Let p : Z ^ MU {-l-oo}, not identically {+00}, be given. We define its subconju- 
gate : X ^ MU {+00} by: 

(2.4) (a;) = sup {u (x, z) — p (z)} 

z 

It follows from the definitions that is a u-convex function on X (it might be 
identically {-l-oo}). 

Example 1. Set f (x) = u {x, z) -f a. Then 

(z) = sup {u (x, z) — u (x, z) — a] — —a 

X 

Conjugation reverses ordering: if /i < /2 , then f\> f\, and if pi < P2 , then p\ > 
p\. As a consequence, if / is u-convex, not identically {-l-oo}, then /" is u-convex, 
not identically {-l-oo},. Indeed, since / is u-convex, we have / {x) > u {x, z) -\- a for 
some (z, a), and then (z) < —a < 00. 

Proposition 1 (the Fenchcl inequality). For any functions f : X MU{-l-oo} 
and p : X —t MU {+00}, not identically {-l-oo}, we have: 

V(x,z), f{x) + f{z)>u{x,z) 
V (x, z) p (z) -f p^ {x) > u {x, z) 
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2.3. Subgradients. Let f : X ^ MU {+00} be given, not identically {+00}. Take 
some point x E X. We shall say that a point z G Z is a subgradient of / at x if the 
points X and z achieve equality in the Fenchel inequality: 

(2.5) f{x) + fHz)^u{x,z) 

The set of subgradients of / at x will be called the subdifferential of / at a; and 
denoted by df {x) . 

Similarly, let p : Z ^ RU {+00} be given, not identically {+00}. Take some 
point z £ Z. We shall say that a point x (z X is a subgradient of p at z if: 

(2.6) (x) + p (z) ~ u (x, z) 

The set of subgradients of p at z will be called the subdifferential of p at z and 
denoted by dp{z). 

Proposition 2. The following are equivalent: 

(1) zedfix) 

(2) Vx', / {x') > f{x)+u {x\ z)-u {x, z) 

If equality holds for some x' , then z G df {x') as well. 

Proof. We begin with proving that the first condition implies the second one. As- 
sume z G df (x). Then, by H2.5I) and the Fenchel inequality, we have: 

/ {x') > u {x, z) - /" (z) = u (a;', z) - [u {x, z) - f {x)] 

We then prove that the second condition implies the first one. Using the inequal- 
ity, we have: 

/« (z) = sup{u(a;',z)-/(x')} 

x' 

< sup {u (x', z) — f {x) — u (x', z) -\- u (x, z)} 

x' 

= u{x,z) ~ f{x) 

so / (x) + f^ (z) < u [x, z). We have the converse by the Fenchel inequality, so 
equality holds. 

Finally, if equality holds for some x' in condition (2), then / [x') — u {x' , z) = 
/ (a:) — u {x, z), so that: 

Vx", / {x") >f{x)-u (x, z)+u (x", z) 

— f (x') — u {x' , z) + u (x", z) 

which implies that z G 9/ (x'). □ 

There is a similar result for functions p : Z ^ RU {+00}, not identically {-I- 00}: 
we have x € dp (z) if and only if 

(2.7) V(x',z'), _p(z') >p(z)-f-u(x,z') -u(x,z) 

2.4. Biconjugates. It follows from the Fenchel inequality that, if p : Z RU {+c)o} 
is not identically {-I-cxd}: 

(2.8) p"* (z) = sup{w(x,z) -p" (x)} < p(z) 
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Example 2. Set p (z) — u {x, z) + b. Then 

p^'^ (z) = sup {u (x, z) — (x)} 

X 

> u {x, z) — p^ [x) 
= u {x, z) + b — p (z) 

This example generalizes to all w-convex functions. Denote by C„ (Z) the set of 
all u-convex functions on Z. 

Proposition 3. For every function p : Z —t RU {+oo}, not identically {+cx)}, we 
have 

(z) = sup {(^ (z) I < p, e C„ {Z) } 

Proof. Denote by p (z) the right-hand side of the above formula. We want to show 
that (z) =p{z) 

Since p'*" < p and p"'* is u-convex, we must have p"'* < p. 

On the other hand, p is u-convex because it is a supremum of u-convex functions. 
So there must be some J5 C X x R such that: 

p (z) — sup {u {x, z) + 5} 

{x.b)eB 

Let (x, 6) G B. Sincep < p, we have u {x, z)+6 < p (z) < p (z). Taking biconjugates, 
as in the preceding example, we get u (x, z) + 6 < p"' (z). Taking the supremums 
over (x, 6) S S, we get the desired result. □ 

Corollary 1. Letp : Z — > MU {+(X)} be a u-convex function, not identically {+oo}. 
Then p — p"' , and the following are equivalent: 

(1) X e dp{z) 

(2) p(z) +pi' (x) ^u{x,z) 

(3) z e 9p« (x) 

Proof. We have p"" < p always by relation (|2.8|) . Since p is u-convex, we have: 

p (z) = sup {u (x, z) -I- 6} 

ix,b)eB 

for some _B C X x R. By proposition O we have: 

sup {u (x, z) + b} < pW (z) 

(2;,fc)e-B 

and so we must have p = p"*. Taking this relation into account, as well as the 
definition of the subgradient, we see that condition (2) is equivalent both to (1) 
and to (2) □ 

Definition 1. We shall say that a function p : Z RU {+cxd} is u-adapted if it is 
not identically {+00} and there is some {x,b) Cz X x R such that: 

\/z ^ Z, p (z) > u (x, z) + b 

It follows from the above that if p is u-adapted, then so are p", p"" and all further 
subconjugates. Note that a u-convex function which is not identically {-l-oo} is it- 
adapted. 

Corollary 2. Let p : Z ^ RU {+cx3} be u-adapted. Then : 
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Proof. If p is u-adapted, then p^ is M-convex and not identically {+00}. The result 
then follows from corollary ^ □ 

2.5. Smoothness. Since u is continuous and XxZ \s compact, the family {u (x, •) \ x £ X} 
is uniformly equicontinuous on Z . It follows from the definition l2 . 2l that all u-convex 
functions on Z are continuous (in particular, they are finite everywhere).. 

Denote by k the upper bound of WD^u (x, z)\\ for {x, z) £ X x Z. Since D^u is 
continuous and X xZ is compact, we have k < 00, and the functions a; — > u (x, z) are 
all fc-Lipschitzian on X. Again, it follows from the definition ?? that all u-convex 
functions on X are /c-Lipschitz (in particular, they are finite everywhere). By a 
theorem of Rademacher, they are differentiable almost everywhere with respect to 
the Lebesgue measure.. 

Let f : X ^ R he convex. Since f — f'^\ we have: 

f [x) = sup {u (x, z) - /» (z)} 

z 

Since /" is w-convex, it is continuous, and the supremum is achieved on the right- 
hand side, at some point z G df {x). This means that all w-convex functions on X 
are subdifferentiable everywhere on X. 

Let X be a point where / is differentiable, with derivative D^fix), and let 
z £ df {x). Consider the function ip {x') = u {x' , z) — /" (z). By proposition [SJ we 
have f < f and (p (x) — f (x) , so that (p and / must have the same derivative at 
X : 

(2.9) DJ{x) = D:,uix,z) 

By assumption (see (jl.HII . this equation defines z uniquely. We shall denote it 
by z = V„/(a;). In other words, at every point x where / is differentiable, the 
subdifferential df {x) reduces to a singleton, namely {V„/ (x)}. Combining all this 
information, we get: 

Proposition 4. For every u-convex function f : X ^ R, there is a map V„/ : 
X —^ Z such that, for almost every x : 

D^f [x) = Dxu{x,z) ^=4>z = V„/(x) 

The following result will also be useful: 

Proposition 5. Let p : Z ^ R be u-adapted, and let x £ X be given. Then there 
is some point z £ dp^ {x) such that p (z) = p*"* (z). 

Proof. Assume otherwise, so that for every z £ dp^ {x) we have p"" (z) < p{z). For 
every z £ dp^ (x), we have x £ dp^"^ (z), so that, by proposition |21 we have 

pW (z') > u {x, z') - u {x, z) + pW (z) 

for all z' £ Z, the inequality being strict if z' ^ dp^ [x) . Set p^ (z') = u (x, z') — 
u {x, z) + p"* (z). We have: 

z' ^ ap« (.t) =^ ip, (z') < pW (z') < p (z') 
z' e dp* (x) =^ p, (z') < pW (z') < p (z') 

so that pz {z') < p (z') for all (z, z'). Since Z is compact, there is some e > such 
that ipz {z') + s < p (z') for all (z, z'). Taking the subconjugate with respect to z'. 
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we get: 

(x) < sup {u {x, z) - ip^ (z')} - £ 

z' 

= sup {u {x, z') — u (x, z') + u (a;, z) — p"' (z)} — e 

z' 

= u (x, z) — p"* (z) — e = (.t) — £ 
which is a contradiction. The result follows □ 
Corollary 3. If x is a point where is differentiable, then: 

(2.10) p(V„p»(x)) =p««(V„p«(a;)) 
and: 

(2.11) pi {x) = u (x, V„p« {x)) ~ p (V„p« {x)) 

Proof. Just apply the preceding proposition, bearing in mind that dp^ {x) contains 
only one point, namely Vup' (x). This yields equation (|2.10() Equation H2.11() follows 
from the definition of the subgradient: 

p« [x) = u (x, V„p« (x)) - (V„p« (x)) 

and equation (|2.1U|) . □ 

2.6. w-concave functions. Let us now consider the duality between Y and Z. 
Given w : y x Z — > i?, we say that a map g :Y ^ RU {— cxi} is v-concave iff there 
exists a non-empty subset A C Z x M such that: 

(2.12) Vyer, 5(y)= inf {t- (y, z) + a} 

and a function p : Z ^ KU {— oo} will be called v-concave iff there exists a non- 
empty subset B C X X M. such that: 

(2.13) P{z)^ inf {viy,z) + b} 

(x,b)eB 

All the results on u-convex functions carry over to u-concave functions, with 
obvious modifications. The superconjugate of a function g : Y RU {— oo}, not 
identically {— oo}, is defined by: 

(2.14) g' {z)^-ini{v{v,z)-g{y)} 

y 

and the superconjugate of a function p : Z ^ RU {— oo}, not identically {— oo}, is 
given by: 

(2.15) p^(y) = inf{t;(y,z)-p(z)} 

z 

The superdifferential dp^ is defined by: 

9p^ iu) = argmin {i; {y, z) — p (z)} 

z 

and we have the Fenchel inequality: 

p{z)+p\y) <v{y,z) y{y,z) 
with equality iff z G dp^ (y) . Note finally that p^^ > p, with equality if p is t;-concave 
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3. The dual optimization problem. 
Denote by A the set of all bounded function on Z: 

p G A <^=^ sup \p {z)\ < 00 

z 

and consider the minimization problem: 



(3.1) inf 



[ pHx)dn- f p\y)dv 
Jx Jy 



Proposition 6. The minimum is attained in problem (P) 
Take a minimizing sequence p„: 

/ pi{x)dii- [ pliy)du^mi{P) 
Jx Jy 

Setting = Pn + a, for some constant a. Then = — a and Qn = Pn~ o-- Since 
/U {X) = v (Y), we have: 

/ qi{x)df,- j qi{y)dv= j pi{x)df,- j pi{y)dv^ini{¥) 

Setting a = — inf^ Pn (z), we find inf^ qn (z) = 0. So there is no loss of generality 
in assuming that: 

Vz, inf Pn (z) = 

z 

which we shall do from now on. 

The sequences pf^ is fc-Lipschitzian. Since p„ > 0, we have 

Pn (x) < sup u {x, z) < max u{x,z) 

Choose Zn such that pn (zn) < 1. Wc then have: 

pi (x) > u {x, Zn) - Pn (^n) > min u {x, z) - I 

x,z 

So the sequence p| is uniformly bounded. By Ascoli's theorem, there is a uni- 
formly convergent subsequence. Similarly, after extracting this first subsequence, 
we extract another one along which converges uniformly. The resulting subse- 
quence will still be denoted by p„, so that: 

Pn^ f uniformly 
pIi^ 9 uniformly 

Taking limits, we get: 

(3.2) / nx)dii^ ( g{y)dv = mi{V) 

Jx Jy 

It is easy to see that pj]* and p^^ ^ 9^ uniformly. Since p^" < Pn and 

Pn >Pn, we have p^^ < Pn, and hence /' < g^ Set: 



P=i(/> + /). 
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bmce P and are Lipschitz continuous, so is p. Since < p < , we must 
have p« < < f and p^ > g^^ > g. Hence: 

pHx)dti~ f p^iy)di^< [ f{x)dii- f g{y)diy = M{P) 
X Jy Jx Jy 

Since ^ P ^ g^ , both sides being continuous functions, the function p must 
be bounded on Z, and the above inequahty shows that it is a minimizer. The proof 
is concluded. 

4. Proof of the main result. 
Let us now express the optimality condition in problem (P). Set: 

s (x) = V„p» (x) 
t{y) = V,p' (y) 

where the gradient maps V„p'* and V„p^ have been defined in proposition^ 
Proposition 7. s(/i) = t (z^) 

Proof. We follow the argument in Carlier ??. Take any continuous function ip : 
Z ^ R. Since p is a minimizer, we have, for any integer n: 



(4.1) 



X 



p+-(p 
n 





c?/i — y n 




) 











diy>0 



We deal with the first integral. Set p + ^ip — Pn- Since pjj is u-convex, it differen- 
tiable almost everywhere. Take a negligible subset Xq such that all the Pn,n € N, 
and p, are differentiable at every x ^ Xq. If a; ^ Xq, then Vupjj (x) is the only 
point in dp^^ (x), and we have, by corollary |21 

u (x, V^pI (x)) - p {y.,pl (x)) < p« (x) = u (x, VuP* (x)) - p {VuP* (x)) 

so that: 

(4.2) u (x, VuPi {x)) - p {Vupi (x)) ~ u (x, V„p« (x)) + p (V„p« (x)) < 
yielding: 

(4.3) pI (x) + (V„p1 (x)) - p« (x) < 
From the definition of p^^, we have, using corollary |31 again: 

u (x, V„p« (x))-p (V„p« (a^))-^'/' (VuP* (x)) < rf, (x) = w (x, V„p«, {x))~p (V„p«, (a;))-^'/5 (V„p1 I 

Rewriting this, we get: 
(4.4) 

{Vupi (x)) - (V„p« (x)) < u (x, VuPi (x)) - p {VupI (x)) - u (x, V„p« (x)) + p {VuP^ {x)) 

(4.5) = pi (x) + {VupI (x)) - p« (x) 

n 

yielding: 

(4.6) ^i</p (V„p« (x)) < (x) - p« (x) 
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Now let n ^ oo. Using corollary 13 we have: 

U {x, VuPI (x)) = Pn {VupI (x)) + pi [X) 

Since Z is compact, the sequence V,jp^ [x) G Z has a cluster point z, and since p„ 
and converge to p and uniformly, we get in the limit: 

u (x, z) — p (z) + (x) 

so that z G dp^ (x). But that subdifFerential consists only of the point V„p* (x), so 
that z = Vup" (x). This shows that the cluster point z is unique, so that the whole 
sequence must converge: 

VupI (x) ^ V„p« (x) 
Inequalities () and () together give: 

(4.7) (V„p« (x)) < n {pi {x) ~ p» {x)) < if (V„p1 {x)) 

Taking limits in the inequalities (|4.3I) and 14.61) . we get: 

Vx ^ Xq, lim n {pi {x) - p« {x)) = (V„p* (x)) 

Similarly, we have: 

Vj/ ^ >o, lim n (pi (x) - p^ {x)j = -tp {VuP^ (a;)) 

where Yq C F is negligible. 

Because of (I4.7|) . we can apply the dominated convergence theorem to inequality 
We get: 

- y (Vp* (x)) dn + J If (yp' (y)^ du>0 

Since the inequality must hold for —cp as well as (p, it is in fact an equality. In other 
words, for any ip : Z ^ R with compact support, we have: 



(4.8) cp (Vp« (x)) dfi = J^ip (Vp' (y)) 



dv 



and this means that s (/i) — t {v), as announced. □ 

Set s (fi) = t [v) = A. This is a positive measure on Z, not necessarily absolutely 
continous with respect to the Lebesgue measure. 
Applying corollary El we have: 

{x) = u{x,s (x)) — p{s (x)) 

p^ {y) = v{y,t{y)) -PiHy)) 

and hence: 

p^x)dfi- p^{y)diy^ [u{x,s{x)) -p{s{x))]dfi- [v {y , t (y)) - p {t (y))] 
X Jy Jx Jy 



X 



u{x,s{x))dfi- J v{y,t{y))dv - 
u{x,s{x))dp,- / v{y,t{y))dv - 



X 



p (s (x)) d^i- j^p [t (y)) dv 
p{z) dX — / p (z) dX 



(4.9) = / u{x,s{x))dn- / v{y,t{y))dv 
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Let (s, t) be a pair of Borelian maps such that s (/i) — t {v). Then, by the Fenchel 
inequahty:: 



ix 



u{x,s{x))dfi- v{y,t{y))diy < [p^x) + p {s (x))] - (v) + P (t (y)) 



X 



dv 



= j p^ [x) dfi 



X 



X 



p{s{x))dfi - / p{t{y))dy 



The last bracket vanishes because s (/i) = t{v). Applying inequality 14. 91 we get: 



u {x, s {x)) dji 



X 



viy,t{y))diy < / u{x,s{x))d^- / v{y,t{y))dv 



X 



This shows that (s, f) is a maximizer,and proves the existence part of theorem^ 
As for uniqueness, assume that there is another maximizer (s', t'). The preceding 
inequality then becomes an equality: 



X 



: {x, s {x)) d^i~ I V {y, t (y)) diy = / [p« (x) +p{s (x))] d^i- [p^ {y) +p{t {y)) 



X 



dv 



which we rewrite as: 

[p^{x)+p{s{x))^u{x,s{x))\dii+ I v{y,t{y))-p^{y)+p{t{y)) 



X 



diy = 



Both integrands are non-negative by the Fenchel inequality. If the sum is zero, 
each integral must vanish, and since the integrands are non-negative, each integrand 
must vanish almost everywhere. This means that: 



s {x) G dp^ {x) a.e. 
t{y)edp^{y) a.e. 

and since dp^ (x) — {s (x)} and dp^ {y) = {t (y)} almost everywhere, the result is 
proved. 

5. Some consequences. 

We shall now investigate some properties of the function p : Z ^ R. 

Recall that we denote A = s{fi) = t{v). It is a positive measure on Z. Its 
support Supp (A) is the complement of the largest open subset C Z such that 
A = on r2. If for instance fi and i> have the property that the measure of any open 
non-empty subset is positive, then: 

Supp (A) ^s{X) =t{Y) 
Proposition 8. p"** = p = p^^ on Supp (A) 

Proof. We have seen that p"" (s (x)) — p (s (a:)) /x-almost everywhere. Since s (fi) — 
A, this means that p'" (z) — p (z) A-almost everywhere. Since p and p**" are contin- 
uous, equality extends to the support of A. Similarly, p = p^^ on Supp (A), and the 
result follows □ 



Let us illustrate this with an example 
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5.1. The linear case. Suppose X,Y,Z are compact subsets of i?", and we want 
to minimize: 

(5.1) / \\x-six)fd^i- [ \\y-tiy)fd,, 

JX JY 

among all maps (s, t) such that s (/i) = t {v). Developing the squares, this amounts 
to minimizing: 



X d^JL — I y dv 

X JY 



s (x)^ dfi — t (y)^ dv 



x's {x) d/i — / y't {y) dv 



Y 



The first bracket is a constant (it does not depend on the choice of s and t) . The 
second bracket vanishes because s [ji) =t{v). We are left wiht the last one. So the 
problem amounts to maximizing: 

x' s {x) dfj. — y't (y) dv 

IX JY 

and it falls within the scope of theorem^by setting u (x, z) = x' z and v (y, z) — y' z. 
Then u-convex functions are convex in the usual sense, u-concave functions are 
concave in the usual sense. By proposition |H1 p is linear on Supp (A), so we may 
take: 

p (z) = tt'z 

for some vector tt € We then get s{x) by maximizing {x — i:)' z over Z. 

Similarly, we get t [y] by minimizing {y — ir) z over Z . Note that this implies that 
s (X) = i (Y) C dZ, the boundary of Z. Let us summarize: 

Proposition 9. There is a single map (s, i) which minimizes the integral ?? among 
all maps (s, t) such that s (/i) — t (v). It is given by: 

(5.2) s (x) = argmax (a; — tt)' z 

z 

(5.3) f (y) = argmin (j/ — tt)' z 

z 

and the actual value of tt E i?" is found by substituting 1^5. and in the 

integral 1^5. 1\ and by minimizing the resulting function of a. 

This first example is degenerate: the dimension of Supp (A) is strictly less than 
the dimension n of X,Y and Z. Let us now go in the opposite direction. 

5.2. The non-degenerate case. Suppose X, F, Z are compact subsets of i?" . Let 
z belong to the interior of Supp (A), so that z — s (x) — i (?/)and suppose there are 
neignbourhoods Q^, ^y, and of x, y and z such that the restrictions s : Qx ^ ^z 
and f : rij, — > il^ are invertible, with continuous inverses a : Qz ^ and r : 

^z ^ 

The function p then satisfies two partial differential equations on fi^ : a second- 
order equation of Monge- Ampere type, and a fourth-order equation of Euler-Lagrange 
type 

5.2.1. A Monge-Ampere equation. .Write the definitions of cr (z) and t (z), for z G 

(5.4) Dzu{a{z),z)=Dzp{z) 

(5.5) Dzv{t{z),z)^Dzp{z) 



AN OPTIMAL MATCHING PROBLEM 
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Inverting the first equation expresses c (z) in terms of Dzp{z). Inverting the 
second one expresses t (z) in terms of Dzp{z). Substituting into the equation 
o T (/i) = V gives a second-order partial differential equation for p. 
Let us give an example. Consider the problem of minimizing the integral: 

1 



(5.6) 



— ■ ||a; — s (a;) II dx + 



X 



\y~t{y)\\ dy 



over all maps (s,t) such that s (/i) — t{v) . Here a > is a given constant. We 



apply theorem^with u {x, z) 



\\x-z\\ andv{y,z) = \\\y~z\\ 



Assume Vtr^ , ily and 0,^ are as above. Equations (|5.4|l and (|5.5|l become: 



DzP{z) = -a(cr(z) - z) 
Dzp{z) = (r(z)-z) 

So (T (z) = z — ^DzP (z) and r (z) = z + DzP (z). The map o r sends /i on 
We must have: 

[det Dza (z)]^'^ det DzT (z) = 1 
and this gives a second-order equation for p : 



(5.7) 



det [/ + Dl^p (z)] = dot 



I--Dl,piz) 
a 



5.2.2. An Euler- Lagrange equation. Assume that the constraints on p are not bind- 
ing on flz- In other words, there is some e > such that, for every h such 
that \h\ < e, and every C°° function tp with compact support in the func- 
tion Ph ^ P + hip is still u-convex and w-concave. 
Recall that p solves the optimization problem: 



inf 

p 



p« (x) df,- p^ (y) di^ 



X 



and this implies that: 



pl{x)df,- pi{y)d,y> pHx)df,- p^{y)dv 

"^Oy "^Oj; "^Oy 

for every h. Expressing the sub- and superconjugates in terms of the sub- and 
superdifferentials, and taking advantage of the fact that pf," = p = we get, for 
every h : 

(5.8) 

u (ah (z) , z) - Ph (z)] d [cr^^ in)] - [v [rh (z) , z) - ph (z)] d [t^'^ {v)] 



(5.9) 



>/ [^.(a(z),z)-p-(z)]d[a-l(^)]- / ^ (r (z) , z) - p (z)] d [r'^ (^)] 



where ah and Th are defined by: 

(5.10) Dzu{<Jhiz),z)^DzPhiz) 

(5.11) (r,, (z) , z) = i^.Ph (z) 
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Note that for ft, = 0, we have ah = cr, Th ~ t and ph — p, while a (/i) = r (/i). 
Letting /i — > 0, we find that p must satisfy an Euler-Lagrange equation of the fourth 
order. 

Let us illustrate this with example (|5.6() . 

A function p : Z ^ R is u-convex iff -D^^p (z) > ~al for every z, where D^^p (z) 
is the Hessian matrix at z, and it is w-concave iff D^^p (z) < I. So there are many 
functions p : Z ^ R which are both w-convex and u-concave: they must satisfy 
—al < D^^p (z) < / everywhere. 

Assume fl^ji^y and flz are as above, and —al < Dl^p{z) < I on fl^. Then the 
integral to be maximized with respect to p — ph is: 

-1 



z - -Dzp{z),z ] -p{z) 



I--DIp{z) 
a 



dz 



[v {z + D,p (z) ,z)^p (z)] [/ - Dl,p (z)] ' dz = 



^\\DM\'-P 
Za 



I--Di 
a 



l\\D.pf-P 



[i-dLp\ 



dz 



So p, which minimizes the last integral, must satisfy the corresponding fourth or- 
der Euler-Lagrange equation We will not write it down explicitly, although relation 
(15.71) . which is valid for p = p, would introduce some simplifications. 
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